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Abstract. If supersymmetry is observed at the LHC its model parameters can be measured at the elec- 
troweak scale. We discuss the expected precision on the parameter determination, including a proper 
treatment of experimental and theoretical errors. Particular attention is paid to degenerate solutions. Us- 
ing the SFitter framework we perform a bottom-up reconstruction of the unified parameters at the high 
scale, including a full error propagation. 

PACS. ll.30.Pb Supersymmetry - 12.60.Jv Supersymmctric models 



1 Introduction 

While supersymmetry was originally not introduced as a 
phenomenological model targeted at definite shortcomings 
of our Standard Model, it has since developed into the 
most attractive ultraviolet completion of the Standard 
Model at and above the TeV scale. The basis of super- 
symmetric theories is a matching of fermionic degree of 
freedom with bosonic degree of freedom, now given the 
particle content of the Standard Model at the electroweak 
scale [TK5JI3] . This symmetry automatically cures the theo- 
retical problem with a perturbatively unstable fundamen- 
tal Higgs mass in the Standard Model, i.e. the hierarchy 
problem. 

Several specific points make the minimal supersym- 
metric extension of the Standard Model (MSSM) an at- 
tractive ultraviolet completion: while electroweak preci- 
sion data favors a light Higgs boson [4, supersymmetry 
provides such a lightest Higgs boson with a mass less 
than about 140 GeV. To implement fundamental symme- 
tries protecting the proton life time and avoiding flavor- 
changing neutral currents into the MSSM wc usually re- 
sort to R parity. This parity in turn forces the light- 
est supersymmetric particle to be stable. Attributed to a 
'WIMP miracle' the typical relic densities predicted for a 
weakly interacting supersymmetric dark matter candidate 
roughly agree with the observed values [SIE]- 

In addition to these attractive weak-scale properties 
supersymmetry offers another, unique, opportunity: it al- 
lows us to predictively extrapolate a perturbative and 
renormalizable gauge theory to arbitrarily high energy 
scales. While we know that the three gauge couplings do 
not unify in the Standard Model, including supersym- 



metric degrees of freedom at the TeV scale can natu- 
rally lead to such a unification [71[Sllllini[IIl[IllI3] ■ This is 
one of the reasons why fundamental and unbroken super- 
symmetry could live above the unification scale Qgut > 
IQis GeV [ni[ni[THl[T71P^pi?l[^P1[^ . One way of break- 
ing supersymmetry would then be a gravity-driven link be- 
tween a hidden sector and the MSSM, where the scales get 
adjusted to accommodate supersymmetry-breaking mass 
parameters around the TeV scale. Using this setup a renor- 
malization group analysis shows that the electroweak sym- 
metry is automatically broken by running one of the two 
Higgs masses squared to negative values slightly below the 
mass scale of the supersymmetric partners [23] . 

In the coming years supersymmetry should be discov- 
ered at the LHC. An e+e~ linear collider (ILC) with a 
center-of-mass energy of 500 GeV extendible to 1 TeV, 
as it is under study, can shed further light on supersym- 
metry. A formidable task will be to determine the funda- 
mental parameters of supersymmetry from experimental 
measurements. Even if we expect supersymmetry to unify 
at a high scale we should not assume such a unification 
and simply fit the high-scale parameters to experimental 
measurements at the electroweak scale. Instead, the ap- 
propriate though technically challenging question should 
be: do the measured weak-scale MSSM parameters evolved 
to high energy scales unify? Or, if not, can we observe 
relics of a unification, like sum rules stable with respect 
to renormalization group running |24j . 

In principle, grand unification of the supersymmetric 
parameters can be observed from data in the gaugino and 
scalar sectors [25l[26ll27] or in the gaugino sector only [28] . 
The aspect we focus on in this study is the determina- 
tion of the central values and errors [2^130) . i.e. a proper 
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measurement of supersymm.etric unification. The bottom- 
up renormalization group analysis are performed for ex- 
pected LHC measurements and its combination with the 
ILC. Because our analysis depends critically on the knowl- 
edge of all errors, we use the well studied parameter set 
SPSla |3T]. It is interesting to note that the result of the 
fit of the electroweak data, adding b-physics observables, 
the anomalous magnetic moment of the muon [32^ and the 
relic density [33], yields a best-fit point not too far from 
SPSla, a further motivation to study in detail such a pa- 
rameter choice [34,35,36 . Note that by adjusting slightly 
the parameters of SPSla (denoted SPSla' or SPAl in 37 ) 
the relic density can be reduced to agree with the mea- 
surement from WMAP [33] without changing the collider 
observables significantly. The relic density is only a side 
product of our analysis, therefore we stick to SPSla for 
which the experimental error estimates are available di- 
rectly. 

Our analysis relies on the SFittcr framework [291130] . 
The same techniques have been successfully applied to 
other questions, like the determination of the Higgs bo- 
son couplings at the LHC [3S]. The theoretical aspects 
of evolving the relevant parameters from the electroweak 
to the GUT scale and vice-versa we discussed in Sec- 
tion [2] The expected measurements at the LHC and ILC 
are listed in Section (3] followed by a summary of the de- 
termination of the MSSM parameters at the electroweak 
scale in Section |3] In Section [5] we develop the method for 
the determination of the unified parameters at the high 
scale and apply it to the respective LHC and LHC-I-ILC 
measurements . 



2 Theoretical aspects 

For the following studies, the MSSM is defined as a variant 
of the phenomenological MSSM with the following param- 
eters (evaluated for convenience at the electroweak scale, 
1 TeV [37]): tan/3 is the ratio of the vacuum expecta- 
tion values of the Higgs doublets. Mi, M2 and M3 are the 
gaugino mass parameters, fi is the higgsino supersymmet- 
ric mass parameter and ttt-a the physical (pole) mass of 
the pseudo scalar Higgs boson. For simplicity the gaug- 
ino masses are restricted to be positive. The soft breaking 
masses in the slepton sector are denoted as Mg^, Mgj^, 
M^l: M^j,, Mf^ and M^^. For the squarks, as u, d, s and 
c quarks are experimentally practically indistinguishable, 
a common definition (average) is used for the first two 
generations denoted as Mq^ and Mq^. The third genera- 
tion is treated separately with the parameters Mqs^^ , Mj 
and Mr . The tri-linear terms, irrelevant for the first two 
generations due to the smallncss of the fermion masses, 
are taken into account for the third generation with A,-, 
At and Af,. 

The supersymmetric and soft-supersymmetry break- 
ing parameters of the MSSM can in principle be defined 
at some arbitrary scale, not only the electroweak scale. In 
order to compare the results of the bottom-up and top- 
down approach, the high scale MSSM will also be used. 
The parameters of this model are defined at the GUT scale 



(■~ 10^^ GeV). In this model the two weak scale parame- 
ters fi and mA are replaced with the soft supersymmetry- 
breaking terms for the Higgs doublets, Mhi and Mhs ■ 



2.1 Renormalization Group Evolution properties 

The renormalization group equations (RGE) [41] play an 
important role in the analysis by relating low to high 
scale MSSM parameters. Corrections up to the level of 
two-loops are implemented in SUSPECT [iU] and Soft- 
SUSY [52]. Three-loop results for the beta functions are 
known and have been shown to stay within the two-loop 
error bands [43ll44] . Unless stated otherwise the full two- 
loop corrections are applied in the following. An approx- 
imate analytical form is useful to understand the RGE 
dependence of the most relevant parameters. Typically, 
at leading one-loop order the RGE for the gaugino mass 
parameters are closely related to the ones for the gauge 
coupling [41] : 
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where i = InQ, Bi = (33/5, 1, —3) andM^, gi fori = 1,2,3 
are the gaugino masses and gauge couplings (in the stan- 
dard normalization with gi — \/5/3gY)- This immedi- 
ately implies 



M,(g) 
gUQ) 



= constant 



^ M,(g,„) 
~ gfiQtn) 



(2) 



where Qin is an arbitrary initial scale, either the high scale 
for top-down or the low scale for bottom-up evolution. 
More explicitly, the one-loop RGE solutions for the gauge 
couplings lead to 



1 + ^adQ^n) In(QVQin') 

47r 



M,(Q) = ZM^{Q^n) 



(3) 



where a^ = gf /(47r). 

The structure of the RGE thus shows that the gaugino 
sector is essentially determined by gauge couplings, and to 
some extent by the Yukawa couplings which only enter at 
the two-loop level. The gaugino parameters are decoupled 
from the scalar sector. In contrast, the soft parameters in 
the scalar sector are strongly correlated due to the RGE: 
most of the soft scalar masses depend not only on their 
value at the initial scale, but also on other scalar masses 
as well as gaugino mass and trilinear coupling parameters. 
An approximation of those solutions for all possible scalar 
masses can be written e.g. in terms of high scale universal 
mQ,mi/2, Aq parameters, but these are not very useful 
and can even be misleading. The RGE for the slepton 
and squark parameters of the first two generations, Mgj^, 
MqL, Mqj^, . . . , depend essentially on the gauge couplings, 
gaugino masses and Tr[yTO^] at one- loop [H]. In fact. 
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type of 
measurement 



nominal stat. LES JES theo. 
value error 



mt 



"^9H - mo 



m-g 



mi- 



's "^2 
max 
HI ■ 
^max. 

mTrixl): 

high 
^Iq ■■ 

, thres 

muq : 

, thres 

mm : 



three-particle edge(x2,^~fl,Xi) 
three-particle edge((7L,X2,xS) 
three-particle edge(gL,X2)^-R) 
three-particle edge(x4,^i,Xi) 
three-particle edge(x2,'?i,Xi) 
four-particle edge(gL,X2,^~fl,X?) 
threshold(gL,X2,^"H.Xi) 
threshold(&i,X2,^~fl,Xi) 



108.7 


0.01 


171.20 


0.01 


102.38 


2.3 


511.38 


2.3 


446.39 


10.0 


89.01 


1.5 


62.93 


2.5 


80.852 


0.042 


449.08 


1.4 


326.32 


1.3 


277.36 


3.3 


83.21 


5.0 


390.18 


1.4 


216.00 


2.3 


198.41 


5.1 



0.25 



0.1 



0.08 



0.3 



1.0 

6.0 
4.3 
1.0 
0.7 

4.3 
3.0 

0.8 
3.8 
2.0 
1.8 



2.0 

1.1 
6.1 
5.5 
8.0 
8.2 
1.2 
5.1 
5.2 
2.0 
1.0 
5.0 
3.3 
3.1 



Table 1. LHC measurements in SPSla, taken from [39]. Shown are the nominal values (from SUSPECT [40]) and statistical 
errors, systematic errors from the lepton (LES) and jet energy scale (JES) and theoretical errors. All values are given in GeV. 
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with Tr[ym^] defined as: 
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i.e. the sum over all scalar soft terms degrees of freedom 
weighted by their hypercharge. This trace has the property 
of vanishing at tree-level for any model where some (even 
partial) universality relations holds among the soft masses, 
and moreover remains constant at one-loop level when 
evolved to an arbitrary scale Q (i.e. ^[Tr[ym^]i_ioop] = 
0)^5]. At two- loop RGE order, it gives roughly a rela- 
tive correction of at most 10% of the largest scalar mass 
squared. In SPSla this squared mass is M^^ . If the Tr[Fm^] 
is zero at one-loop level, i.e., all scalar parameters are well 
determined in SPSla, a moderate increase of the errors 
on the parameters is expected after their evolution to the 
high scale. Due to the relatively large coefficient of the 
asMi term in Eq. [5] the RGEs are sensitive to M3 . 

For the third generation scalar masses, the RGEs are 
more involved and definitely couple different scalar species. 
Typically for the relevant parameters Mf^ , ^fu. j their re- 
spective (one-loop) RGE both involve (apart from Tr[ym^]) 
the parameters: {Mhi, Mf^, Mf^, At-}, plus the relevant 
gauge couplings and gaugino masses. Thus for instance 
the parameter Mf j^ of the MSSM at the final (low or high) 
scale after RGE running from the initial scale will depend 
on the value of the initial Mf^ value, as well as the other 
parameters above, in a complicated way. 



To illustrate the impact of these correlations, the pre- 
cision of the determination of the MSSM parameters will 
be compared to the precision of the high scale MSSM pa- 
rameters. 



3 Collider Data 

The SPSla parameter set leads to moderately heavy squarks 
and gluinos in the range of 500 GeV to 600 GeV. The slep- 
tons have masses ranging from 130 GeV to about 200 GeV. 
The light neutralinos and chargino have masses well below 
200 GeV and their field content is predominantly gaugino, 
whereas the heavier states are higgsino. The detailed anal- 
yses at the LHC and the ILC can be found in ReL [39] . 



3.1 LHC and ILC measurements 

At the LHC the SPSla spectrum can lead to long decay 
chains, the most prominent being: 



QL ^ x°2q -^ ^R^q ^ iiqx'i- 



(7) 



In this decay chain at least five observables can be deter- 
mined P51l47| . The observables are endpoints or thresh- 
olds of invariant mass combinations among the leptons 
and jets. Additional measurements cover essentially the 
strongly interacting sector. In SPSla it is difficult to ob- 
serve the stop quarks above the supersymmetric back- 
ground from sbottom decays leading to the same final 
state. The use of stop sector branching ratios has been 
explored in [551H5] . Stop kinematic edges have been stud- 
ied for other benchmark points [^[50] . While the corre- 
sponding results have not been experimentally confirmed, 
recent progress in fat-jet analysis techniques indicates that 
by the time LHC acquires a sufficient luminosity we should 
be able to measure the stop mass as well |51j . 
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Parts of the electroweak sector, namely three of the 
four neutralinos, Xi, X2 ^-^d xl but not X3 will be observed 
at the LHC. The absence of the fourth neutralino leads 
to ambiguities, e.g., one could suppose that Xs ^as been 
measured instead of xi- In Ref. [3^ an example of the 
consequences of such a wrong assignment is discussed. In 
the following such discrete ambiguities will be left out of 
the discussion. 

In the slepton sector the first and second generation 
selectrons and smuons will be measured as well as the 
lightest stau. The expected precision at the LHC for the 
measurements is listed in Table [T] for an integrated lumi- 
nosity of 300 fb~^. In the analysis, the channels involving 
leptons have been separated for muons and electrons. The 
systematic error of each channel was unchanged, but the 
statistical error was increased to take into account the re- 
duced statistics per observation. This approach should be 
considered as the optimistic limit of what can be done at 
the LHC as additional systematics, e.g., due to a more 
difficult fit of the background, would have to be added. 
Note that even with the increase of the statistical error, 
the systematic (energy scale) error still dominates the ex- 
perimental error. Additional observables such as the use 
of cross sections times branching ratios have been studied 
in [52] • The masses of the sparticles can be derived from 
the observables listed in Table [T] with a fit or toy experi- 
ments without the use of the underlying theory ,39, . 

At the ILC essentially all kinematically accessible spar- 
ticles can be measured. Masses are measured either in 
direct reconstruction at a center-of-mass energy higher 
than the production threshold or via a measurement of 
the production cross section as function of the center-of- 
mass energy at the threshold of (s)particle production. As 
the beam energy is well known from the accelerator, typ- 
ically the expected precision of the mass measurements is 
about an order of magnitude better than at the LHC. The 
ILC observables are shown in Table [2] 

As the RGE running depends strongly on the top quark 
Yukawa coupling value and its error, the pole mass of 
the top quark and the strong coupling constant as are 
included as parameters and measurements in the fit, in 
addition to the supersymmetric and soft-supersymmetry- 
breaking parameters. An error of 1 GeV is assumed when 
only LHC data is used [531HU] . For parameter determina- 
tions involving the ILC, a theoretical error of 0.12 GeV is 
used with a negligible statistical error [SJ • For the strong 
coupling constant a conservative error estimate of 0.001 is 
used [55j. 



3.2 Theoretical analysis of the neutralino sector 
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2.0 
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0.7 


Tl 
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Table 2. Errors for the mass determination in SPSla, taken 
from 39 . Shown are the nominal parameter values (from SUS- 
PECT 40 ), fixing the electroweak symmetry breaking and 
renormalization scales at 1 TeV, the error for the ILC alone as 
well as the theoretical error, all in units of GeV. 
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(8) 



where sw = sinOw, cw = coa0w, s/3 — sin/3, c^ = cos/3. 
The dominant radiative corrections to the neutralino masses 
[57] are incorporated in the form of Mi -^ Mi + Z\Mi, 
M2 ^ M2 -I- Z\M2, fi ^ n + Afi (tan/3, mz and sw 
can als o be considered as the radiatively corrected values 
in the DR scheme). Assuming that the four neutralino 
masses are measured, fixing tan/3 temporarily, the first 
approximation is mz — )■ 0: in this case, the correspon- 
dence between the mass eigenvalues and basic parameters 
of Eq. [5] is trivially given as 



,o(i = l,..,4)=Mi,M2,|At|,|M| 



(9) 



with all possible permutations, i.e. one obtains a 12-fold 
degeneracy, corresponding to the six possible permuta- 
tions among the parameters Mi, M2, |/z| and the ambi- 
guity of the sign of /x. Restoring the full mz dependence 
renders the solution more complex but qualitatively sim- 
ilar: given that the three neutralinos are measured at the 
LHC, Ml, M2 and fj, are determined from the following 



The measurements of the neutralino masses strongly in- 
fluence the determination of the parameters. Most of the 

qualitative results of the full analysis in the gaugino/higgsino ^^^^^^ ^^ ^^^^^ equations m^ 
sector parameters can be understood from a simplified 

theoretical analysis, which depends solely on the neutralino „, ,9 9 t,^^^ ,-.^ 

■^ ' ^ -^ "2 , (^^2 + rn%- M1M2 + (M 



sector parameters, neglecting all errors. The neutralino 
mass matrix at tree level is : 



p2 



M2)S,, ^ Sf^)B,, 



+fimlicl,Mi + sl,M2) sin2/3 - fi^MiM2 = (10) 
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and 



(Ml + M2 - 5,,)4- + {pHMi + M2) 



-771 
,2 f^2 
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M2-Msin2/3))Py 



'^{mlic^Mi 



52,,M2)sin2/3 
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^MiM2)5,, 
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with & 



rn^o + m^o , Pj^ 



m^qm^o for any pair of 



neutralinos i,j — 1,2,4(3, and 

Ai' = MiM2-to|-(Pi24 + 5i24(Mi + M2 - 5i24)) (12) 

where 6*124 = ^fn^o + m^a + to-^^o and P124 = m^om^o + 

Al A2 -^^4 -^.1 A2 

Al X4 X2 A4 

For SPSla this system can easily be solved numeri- 
cally to obtain the full set of degenerate solutions for Mi , 
M2, fi, labeled DSl to DS12 in Table El These solutions 
cover all possible hierarchies among |Mi|, IM2I and |/x|. 
Due to mz ^ 0, they no longer correspond to simple 
permutations. The six different hierarchies remain clear 
and the solutions for the opposite sign of /x are not ex- 
actly symmetrical. In fact the system Eq. [TUl Eq. [TT] and 
Eq-im gives 12 solutions, not necessarily all real- valued, 
for any fixed neutralino mass input. Taking different pos- 
sible (physically irrelevant) sign choices for the input neu- 
tralino masses exhausts all possible solutions of different 
/z. Ml signs within the six different hierarchies. In the 
present study only solutions with positive Mi are consid- 
ered for simplicity. Simple approximate solutions are de- 
rived in Ref. [37] by expanding Eqs. lTUlfT^ to first order in 
m%. The difference with respect to the values in Table [3] 
is about one percent. 
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Ml 


M2 


M 


'^prcd 


DSl 


97.66 


187.35 


-358.43 


367.7 


DS2 


182.44 


98.54 


-361.64 


371.8 


DS3 


102.35 


354.88 


-184.61 


195.5 


DS4 


368.7 


120.16 


-165.49 


197.0 


DS5 


168.0 


357.44 


-115.27 


127.3 


DS6 


369.45 


144.05 


-77.94 


55.2 


DS7 


100.41 


196.68 


349.34 


355.6 


DS8 


184.73 


106.47 


354.69 


361.5 


DS9 


109.26 


350.25 


185.84 


193.2 


DSIO 


367.13 


140.59 


170.86 


215.9 


DSll 


163.59 


354.52 


126.55 


134.6 


DS12 


368.88 


136.13 


83.44 


46.7 



Table 3. The 12 degenerate solutions found in the theoretical 
analysis of the neutralino sector. All values are in GeV. 



All 12 degenerate solutions are compatible with a con- 
sistent radiative electroweak symmetry breaking |/i| solu- 
tion, provided that the values of the other parameters in 
the Higgs sector are calculated consistently from this value 
of fi. 



^ These equations are symmetrical under all neutralino mass 
permutations. 



The values of the remaining neutralino mass (Xmod)' 
uniquely predicted for any of the 12 solutions, are also 
given in Table [31 As expected, Xprcd strictly speaking dis- 
tinguishes the 12 solutions, but it is not measured at the 
LHC. In eight of the solutions the Xpred i^ almost as heavy 
as the X4: but not degenerate in mass due to mz ^ 0. In 
the other four solutions DS5, DS6, DSll and DS12 this 
neutralino becomes the LSP or next-to-lightest supersym- 
metric particle. 



3.3 Errors 

The measurement of unification in the supersymmetric 
sector relies not only on a precise estimation of the ex- 
perimental error, but also on a rigorous treatment of the 
theoretical error. The theoretical error is 0.5% for the 
masses of colorless particles, the neutralinos, charginos 
and sleptons. The error on the gluino and squark mass 
predictions is taken to be 1%. The errors reflect the differ- 
ence between spectrum generators calculating the observ- 
ables with the same precision but using different methods, 
as well as (to some extent) the renormalization scale de- 
pendence as a measure of not yet calculated higher order 
terms. For SPSla, performing with SuSpect a rather large 
variation of the renormalization scale, from 200 GeV to 
1 TeV, gives variations of the Higgs and sparticle masses 
which are comfortably below the quoted uncertainties in 
Table[21 To illustrate that also the difference between spec- 
trum generators is covered will be illustrated by using 
SoftSUSY instead of SuSpect. At the LHC the masses are 
not measured directly, the theory errors are considered to 
be uncorrelated and propagated to the observables. For 
the observables of SPSla this is a conservative choice as 
positively correlated theory errors on masses will lead to 
smaller errors by a factor 2-3. 

The expected precision of the measurements is shown 
in Table [U The last column in this table is different from 
the one shown in Ref. [3^ as the theoretical errors of the 
MSSM predictions are shown here whereas in Ref. [30] the 
errors on the mSUGRA predictions are shown. 

The experimental systematic errors at the LHC such as 
the lepton energy scale are considered to be 99% correlated 
to assure that the correlation matrix can be inverted even 
for negligible statistical errors. For the treatment of the 
theoretical error, the RFit scheme [59] is largely followed. 
Given a set of measurements d and a general correlation 
matrix C 




\d,. 



di\ < a. 
di\ > a. 



(thoo) 

i 

(thco) 



(13) 



where di is the i-th data point predicted by the model 
parameters and di the measurement. The contribution to 
the x^ of a given measurement is zero within one unit 
of the theoretical error. This ensures that no particular 



Claire Adam et al.: Measuring Unification 



parameter value is privileged within the theoretical error 
range. This type of behavior is appropriate for theoretical 
errors as, e.g., higher order corrections will necessarily lead 
to a shift of the prediction within the region of the the- 
oretical error. The shift, contrary to experimental errors, 
has no reason to be distributed like a Gaussian. Outside 
of the theoretical range the experimental error is used. 



4 Parameter determination at the 
electroweak scale 

To find the true parameter set from a set of measurements 
which have statistical errors, theoretical errors and corre- 
lations in a highly correlated system, the SFitter frame- 
work has been developed. SFitter provides several algo- 
rithms to search for the x^ minima (or log-likelihood max- 
ima): weighted Markov chains [5U', a Grid approach and 
a gradient fit (MINUIT). 

SoftSUSY gU and SUSPECT gD] provide predictions 
for the mass spectrum corresponding to the chosen MSSM 
parameters. The mass spectrum is turned into a set of 
expected measurements, which are used as input to the 
MSSM parameter fit. The prediction of the relic density 
is calculated with micrOMEGAS [601I6T] . Unless stated 
otherwise, SUSPECT is used in the following. 

Using the MINOS algorithm from the MINUIT suite 
of tools, the parameters can be determined in a single fit 
together with their errors. However, the analytical propa- 
gation of the errors as function of the scale is quite tedious 
and not always possible, therefore toy experiments are 
used. Typically 5000 (toy) datasets are generated, where 
the expected measurements are smeared according to their 
experimental and theoretical error, including correlations. 
The determination of the parameters is performed for each 
one of the datasets. SoftSUSY and SUSPECT provide 
the RGE running of the fitted parameters between the 
EW and the GUT scale. At any given scale, the width of 
the parameter distributions, either the RMS (Root-Mean- 
Square) or the sigma of a Gaussian fit, is the error on the 
parameter, the mean is the central value of the parameter. 
From here on a parameter set is defined to be the best-fit 
result for a given toy dataset and the RMS is quoted as 
error. 



4.1 Parameter determination from LHC observables 

The number of observables at the LHC is smaller than 
the number of supersymmetric parameters to be deter- 
mined. Therefore two parameters, for which the LHC has 
small or no sensitivity, are fixed. Fixing the parameters to 
the true SPSla values is a solution which can be justified 
a posteriori when grand unification has been proved. In 
this study we have taken a more conservative approach 
of deliberately fixing two trilinear couplings. A,- and A^ 
to the central value of the allowed parameter range, i.e., 
to GeV. Of course, by fixing parameters in a correlated 
system, the errors on other parameters are reduced artifi- 
cially. 



The fixing of the two parameters is de facto a shift of 
250 GeV in A,- and of 750 GeV in A;,. Using the nominal 
values of the other parameters and the true dataset with- 
out smearing but with theory errors, the x^ remains at 
zero. However, performing the same calculation without 
the theory errors, the x^ is 0.8. Thus indeed the two fixed 
parameters have only a small impact on the prediction of 
the LHC observables. The two major contributions to the 
X^ are the edges involving the sbottom masses and the 
gluino mass. As M3 as well as the squark mass parame- 
ters of the third generation are free, these can compensate 
the shift in the prediction introduced by fixing Af,. Since 
the sbottom masses are also used in measurements which 
involve neutralinos and right-handed sleptons of the first 
generation, other parameters such as tan/3. Mi and M2 
are also affected. The magnitude of the shift will be dis- 
cussed later. Suffice it to say at this stage that the effective 
shift of the parameters depends also on the effective corre- 
lation of the measurements, either via the explicit correla- 
tions (energy scale) or that introduced by the flat theory 
errors. The central values of the fits with and without the- 
ory errors are therefore not expected to be the same. The 
latter expectation is verified by defining the theory errors 
as Gaussian systematic errors and compared to the case 
where no theory errors are used. The results for Mi and 
M2 differed by about 1 GeV in the two cases. 

Using Markov chains in the analysis of the weak-scale 
parameters of the MSSM, as pointed out in Ref. |30) . 
eight degenerate solutions are observed with the LHC data 
set. They cannot be distinguished from each other via 
the analysis of the x^ of the best-fit result as they are 
zero when theoretical errors are included. While the hig- 
gsino/gaugino sector is violently different and distinct, 
the other parameters are shifted only slightly between the 
eight solutions. 



The characteristics described in section 13.21 are ob- 
served : for each solution with a positive sign of the hig- 
gsino mass parameter /i, there is an approximate mirror 
solution for negative /z, as well as the permutations of Mi, 
M2 and |^|. 

Note the absence of the four expected solutions DS5, 
DS6, DSll and DS12 corresponding to the "higgsino LSP" 
hierarchies {\^\ < (Mi,M2)) which drastically change the 
neutralino mass hierarchies. A mass splitting of more than 
about 40 GeV between the lightest two neutralinos cannot 
be achieved in this scenario, but a mass splitting twice as 
large is necessary for the LHC observables: the typical 
value of the x^ is of the order of 10^. 

In Tabled the gaugino-higgsino sector is shown for all 
eight degenerate solutions. The numbering defined in Ta- 
ble [3] has been kept to allow comparisons, DS7 is SPSla, 
i.e., the true solution. The errors on the parameters are 
within 20%. These degenerate solutions are indeed well 
defined local minima from which a simple gradient fit like 
MINUIT cannot escape. The central values in Table U 
agree well with the theoretical analysis, which neglects 
all errors, typically to better than one standard deviation, 
validating the results of the two analyses. 
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DSl 


DS2 


DS3 


DS4 


DS7 


DS8 


DS9 


DSIO 


tan/3 


12.3±5.6 


12.4±5.0 


14.9±9.8 


8.9±5.9 


13.8±7.5 


12.6±7.9 


19.2±14.3 


23.0±15.6 


Ml 


102.7±7.1 


189.5±6.2 


107.2±9.2 


383.2±9.1 


105.0±6.9 


191.7±6.6 


116.3±7.5 


380.9±9.3 


Ma 


185.5±7.0 


96.±6.4 


356.9±8.7 


114.2±10.7 


194.7±7.3 


105.5±7.3 


354.0±8.2 


137.2±9.1 


M 


-362.7±7.8 


-364.7±6.8 


-186.0±8.5 


-167.0±9.6 


353.0±7.7 


357.1±8.3 


188.9±7.1 


172.8±8.7 


^xLc 
ILC 


73 

n 

0.17 ±0.07 


22000 

xf 

(4±2) ■ 10"'' 


1700 

xi 

0.14 ±0.08 


25000 

xf 

(8 ±4) -10-* 


0.4 

n 

0.16 ±0.07 


22000 

xf 

(4 ± 3) • lO--* 


2000 

X§ 
0.11 ±0.06 


24000 

X? 
(9 ±4) •10"'' 



Table 4. The result of the parameter determination in the gaugino-higgsino sector is shown for the eiglit degenerate solutions at 
the LHC, including theory errors. DS7 is the true solution (SPSla). The increase of the x^ when adding the ILC measurements 
is shown together with the dominant source of the increase. The last line is the i7h^ prediction from the LHC measurements. 
AH masses are in GeV. 



4.2 Parameter determination from LHC+ILC 
observables 

The addition of the ILC measurements allows to lift the 
degeneracy of the LHC ambiguous solutions easily. Table|4] 
shows the increase of the x^ due to the ILC. Theoretical 
errors are included and no smearing is performed. 

The second-to-last line of Table S] shows the ILC mea- 
surement with the largest contribution to the x^. The 
inversion of Mi with M2 in DS8 (and DS2) is excluded 
via the chargino mass measurement which is sensitive to 
the value of M2. In DS9 (and DS3), where /i and M2 
are exchanged, the chargino mass measurement is not the 
most sensitive measurement as the chargino mass matrix 
is blind to the interchange of these two parameters. Only 
the deviation of their values from an exact exchange leads 
to some sensitivity (10% of the log-likelihood increase). 
Here the third heaviest neutralino, not measured at the 
LHC, leads to a clear distinction with respect to the true 
solution. For DSIO (as well as DS4) again the chargino 
mass measurement, sensitive to the values of M2 and /i 
provides the most stringent distinction from DS7. 

The log-likelihood is calculated using the LHC param- 
eter set where the trilinear couplings A(, and A^ are fixed 
to zero, therefore an increase of the log-likelihood is also 
observed for the true parameter set (DS7) where the effect 
of increasing A,- from —250 GeV to zero is observed via the 
mixing in the precisely measured fi mass. In DSl, which, 
to first order, differs from the true set only in the sign of 
/i, the Ti mass measurement via the mixing provides the 
highest sensitivity. 

The results of the determination of the parameters are 
shown in Table [H] for the LHC and for the LHC combined 
with the ILC. As discussed in the beginning of the section, 
for the LHC the gaugino masses are shifted slightly by 1- 
2 GeV with respect to the nominal value to compensate 
for the fixing of A;, and A,- . Parameters with large errors 
also contribute to this shift. However these shifts are small 
compared to the errors of typically 7 GeV. 

The difference of the results listed in Table [S] with re- 
spect to the previous publication are the following: for the 
LHC the MSSM errors discussed in Section [3] are used in- 
stead of the mSUGRA errors. Additionally the degenerate 
solutions for Aj shown in Ref. ^30] are not separated out 





LHC 


i 


LHC±ILC 


SPSla 


tan/3 


13.8± 


7.4 


10.7± 


3.1 


10.0 


Ml 


105.0± 


6.9 


103.1± 


0.7 


103.1 


Ma 


194.7± 


7.3 


193.0± 


1.6 


192.9 


Ms 


568.3± 


11.6 


568.5± 


7.8 


567.7 


Mf, 


321.4± 


248 


192.4± 


4.7 


193.5 


M.R 


164.3± 


120 


134.9± 


5.7 


133.4 


Ma. 


196.3± 


7.6 


194.4± 


1.2 


194.3 


Mak 


138.0± 


7.0 


135.8± 


0.6 


135.8 


Me. 


196.4± 


7.5 


194.3± 


0.8 


194.3 


Me, 


137.9± 


7.1 


135.8± 


0.6 


135.8 


Mq3. 


491.4± 


16.2 


486. 2± 11.1 


481.1 


Mt^ 


483.4± 


232 


409.6± 17.1 


409.4 


Mi;, 


502.6± 


15.3 


499.1± 13.1 


502.7 


Mq. 


529.6± 


12.1 


526.4± 


5.3 


526.4 


M^a 


508.9± 


16.4 


507.8± 14.4 


506.8 


A, 


fixed 


-102.9± 


681 


-249.3 


At 


-394.4± 


353 


-497.3± 


74 


-496.8 


Ai, 


fixed 


-274.2±1830 


-764.0 


niA 


558.2±271.2 


394.9± 


1.5 


394.9 


^^ 


353.1± 


7.7 


350.8± 


2.5 


351.0 



Table 5. Results for the general MSSM parameter determina- 
tion in SPSla using flat theory errors. The kinematic endpoint 
measurements are used for the LHC and the mass measure- 
ments for the ILC. The LHC+ILC column is the combination 
of the two measurements sets. Shown are the nominal param- 
eter values and the result after fits to the different data sets. 
The MSSM theory errors are used. All masses are in GeV. 



leading to a larger error on At and tan/3. For the combi- 
nation of LHC and ILC, the Higgs mass measurement of 
the ILC is used instead of the LHC measurement. These 
changes are reflected in the significantly smaller errors on 
the parameters. 



4.3 Relic density and Tr[Ym?] 

Any observable sensitive to the neutralino couplings 
and its actual Wino, Bino and Higgsino content, rather 
than only the mass, can help to disentangle the LHC de- 
generate solutions. The relic density fih? for a neutralino 
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LSP is extremely sensitive and drastically changes e.g. for 
(Ml, M2) exchanged hierarchies. But ilh'^ is less sensi- 
tive to the (|/^|,Mi) exchange and even less to the fi sign. 
While the detailed analysis of the relic density is beyond 
the scope of this paper, note that two distinct populations 
can be identified among the eight degenerate solutions. 
The SPSla relic density of 0.19, a factor 1.7 too high with 
respect to the WMAP measurement of 0.1109±0.0056 [33,, 
is obtained for the solutions where the Mi is the smallest 
parameter, i.e., in DSl, DS3, DS7 and DS9. The lightest 
neutralino is essentially a Bino in the bulk. In all other 
cases the relic density is off by three orders of magnitude. 

It is interesting to note that the trace expression de- 
fined in Eq. [5] plays an important role in the stabilization 
of the results. As discussed in section I^TTI this trace is zero 
up to one-loop in models with universality in the scalar 
sector, rendering the first two generation sfermion masses 
very mildly dependent on other scalar terms than them- 
selves. The two-loop corrections lead to a non zero value 
(~ —1.3-10^), which in relative units is a moderate pertur- 
bation within the evolution of most of the scalar masses, 
except for e.g. Mg^^ where it is a substantial contribution 
to its RGE, see Eq.m To test the impact of this additional 
constraint this trace is required to be compatible with its 
SPSla value within 10%. Technically Ti[Ym'^] is added to 
the LHC observables as an additional observable with a 
Gaussian error of 10%. The RMS of the stau parameters 
and the uia is reduced by a factor 5 to 10. 

A consequence of Tr[yTO^] ~ is that in the poorly 
determined third generation the requirement on Tr[yTO^] 
will prevent large values of the stau sector parameters, 
thus reduce the allowed space for these parameters. The 
strong reduction of the error on the poorly determined 
parameters shows the sensitivity of this single constraint. 
From here on the constraint is not used in any of the 
studies. 

5 Evolution to the high scale 

The determination of unification of the supersymmetric 
parameters for the true central values of SPSla a priori 
does not need any special treatment. However in a real 
experiment the measured parameters will be shifted from 
the true values within their error. Due to the coupling in- 
troduced in the RGE equations, some badly measured pa- 
rameters will strongly affect the convergence, in particular 
at the LHC. Therefore there are two separate questions to 
be answered which are intimately related. The first one 
is whether there is a unification of the N parameters and 
second question is what is the value of the unified param- 
eter and at which scale is the unification observed. 



5.1 Bottom-up evolution 

The evolution from the low scale to the high scale is per- 
formed in the following way: for each toy dataset the 
best-fit parameter set is determined at the EW scale. The 
range from 1 TeV (where the parameters are defined) to 
3 • 10^^ GeV (beyond the scale where grand unification is 
expected) is divided into 1000 logarithmically equidistant 



steps. Using SUSPECT [30] the parameters are evolved 
from the EW scale to the next scale point. For each of 
the toy experiments, the fundamental parameters are then 
known at 1000 discrete scale points. 

As far as the RGE are concerned, the evolution of the 
parameters between two scales, for one fixed point in the 
input parameter space, should be independent of whether 
the evolution is performed in a top-down or bottom-up 
manner (apart from negligible numerical integration er- 
rors), as the RGE is obviously invertible. 
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Fig. 1. (left) Mej^ as function of the logarithm of the scale for 
LHC measurements, (right) The fraction of valid non tachyonic 
parameter sets is shown as function of the scale for the LHC 
and the combination of LHC with the ILC. 



However, the errors are amplified strongly as function 
of the scale, especially in the scalar sector, at least for 
some parameters. This is to a large extent a manifestation 
of the "focus-point" phenomena in the MSSM [S2]: even 
if SPSla does not correspond to what is usually referred 
as 'focus-point' scenario in MSSM (which rather corre- 
sponds to much larger toq values), the focusing behavior 
is more general, i.e. in a large part of the MSSM parame- 
ter space the final (low scale) values of some of the scalar 
parameters (in particular Mh2 driving the radiative elec- 
troweak symmetry breaking) are not very sensitive to the 
initial (high scale) input choice. This means that for the 
bottom-up direction even relatively small errors in some 
of the low scale parameters can result in large errors when 
evolved at high scale. Typically it was shown in Ref. [37] 
(Table X in Appendix B) that for relative uncertainties at 
the low electroweak scale of only 1% in the gluino mass M3 
or the up-Higgs doublet mass term Mh2 (letting all other 
parameters at their central SPSla value), the RGE evo- 
lution up to the GUT scale amplifies the errors, resulting 
in relative uncertainties of 20-30% or even 100% on some 
of the final high scale soft mass parameters. The param- 
eters entering Tr[ym^] are those particularly sensitive to 
this divergence behavior. Therefore, if the initial error is 
in the few percent range, some of the sfermion masses can 
become tachyonic well before reaching the final high scale. 
As an illustration Figure [1] (left) shows Msj, as function 
of the logarithm of the scale for all parameter sets. Mg^^ 
is particularly sensitive to the value of Tr[yTO^], as de- 
duced from Eq.[l| Tr[y77i^] can deviate substantially from 
its nominal SPSla value, e.g. from the largely undeter- 
mined M^L in Table [5]) for the LHC, and thus drive it to 
a tachyonic value well before the high scale is reached. A 
strong non-linear scalar mass dependence enters the RGEs 
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of other scalar masses in addition to the Tr[F?7i^], such 
that some tachyonic masses may infect other scalar mass 
RGE. All sets which have at least one tachyonic parameter 
have to be removed. This necessity is also confirmed by 
the analysis of the covariance matrix of the parameters as 
function of the scale. If these tachyonic parameter sets are 
not removed, the covariance matrix can become singular. 

Figure [1] (right) shows the percentage of valid non- 
tachyonic parameter sets. While at the LHC alone the per- 
centage decreases immediately after the electroweak scale, 
the addition of the ILC stabilizes impressively the validity 
of the sets. For the LHC at lO^^ GeV only 30% of the pa- 
rameter sets are still valid, whereas with the addition of 
the ILC 90% remain. At the unification scale for the LHC 
only 7% of the parameter sets remain, whereas for the 
LHC plus ILC measurements 38% remain valid, marking 
a clear improvement over the LHC alone. Similar results 
are obtained using SoftSUSY. 

Once a real measurement is available, toy experiments 
will be defined around the central value of the measured 
data. In the following all confidence level definitions are 
defined with respect to valid, non-tachyonic parameter 
sets. 

Given N parameters for which the grand unification 
is to be tested, the following Xavq i^ to be minimized for 
every scale: 



2 f^2 



N 

(Q^)=V(M, 



mu){Cp%{M,-mu) (14) 



«j 



where Cp is the covariance matrix of the parameters and 
Mj the j'-th mass parameter. 

The scale where the xivn is minimal is the best-fit 
unification scale Qu and the parameter mu is the value 
of the unified parameter. As this procedure is applied to 
each dataset, the resulting distribution of all mjj and Qu 
allow to read off the unification scale and unified param- 
eter value as central values of their distributions and the 
error as RMS or Gaussian sigma of the distributions. 

A closed formula can be derived for the parameter 
mu [63] : 



^um = J2(^p') 



J2iC;%M, (15) 



However this is not sufficient to claim grand unifica- 
tion as these calculations can also be performed for non- 
unifying parameters. To quantify the unification, the ab- 
solute value of Xavg is used. The value is large when the 
N parameters are not compatible with a unified one. It is 
small if the parameters are compatible. If xivg is smaller 
than a cut-off value (xgs): the dataset is unified. The cut- 
off is defined so that 95% of truly unifying datasets have a 
xivg value smaller than X95 • As an example the xivg fo^ ^ 
sample of datasets is shown in Figure [2] as function of the 
scale (DS7). The minimum at a scale of about 10^^ GeV 
is clearly visible. 
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tion of the scale for DS7. The minimum of the distribution is 
observed for a scale of about 10^'' GeV as expected. 



5.2 Evolution of the parameters from LHC observables 

In Table[5]the result of the determination of the MSSM pa- 
rameters is shown in the first column for the LHC. Starting 
from these values the parameter sets are evolved individ- 
ually to the high scale. 

As noted before, at the LHC an eight-fold ambiguity 
will be observed. Therefore the first question is whether 
the RGE evolution of the eight solutions will result in 
similar or different patterns. 
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Fig. 3. Evolution of the gaugino mass parameters to the GUT 
scale for DS7 (SPSla). 



The evolution of the gaugino mass parameters is shown 
in Figure H for DSl, DS2, DS3 and DSIO and in Figure El 
for DS7 (SPSla). The solutions DS8, DS9 and DS4 show 
the same pattern as DS2, DS3 and DSIO, as expected, 
as only the sign of /i changes. In DS2 Mi and M2 are 
exchanged with respect to the correct solution. This leads 
to an intersection of Mi and M3 at 10^^ GeV. In DS9 M3 
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Fig. 4. Evolution of the gaugino mass parameters to the GUT scale for the ambiguous DSl, DS2, DS3 and DSIO at the LHC. 



intersects with M2 at about 10^ GeV, whereas in DSIO 
M3 and Ml intersect at 10^ GeV. 

As expected the correct solution unifies the high scale. 
DSl, qualitatively at least, might unify. Thus six of the 
eight ambiguous solutions can be eliminated as candidates 
for unification. The difference between DSl and the true 
solution being only the sign of /i, it is natural that this 
solution is harder to distinguish from the correct one. 

A comparison of the number of parameter sets com- 
patible with a unified gaugino mass parameter of DSl and 
DS7 is therefore necessary to quantify how well one will 
be able to distinguish the (non)-unification of these two 
solutions. For the true solution (DS7), at the unification 
scale determined for the gauginos, 95.4% of the toy ex- 
periments unify. In DSl only 38% are unified. Thus the 
exclusion of unification for DSl will indeed be very diffi- 
cult at the LHC. 

The results are in agreement with the expectation from 
the structure of the RGBs in the gaugino sector. The ab- 
solute value of the gaugino mass measured at the elec- 
troweak scale gives the starting point of the evolution, 
but the slope is essentially independent of the absolute 
value, so that a wrong parameter value at the electroweak 
scale cannot be compensated. 

As two of the trilinear couplings are fixed at the LHC, 
no further information on the unification can be obtained 
from these parameters. The same is true for the third gen- 
eration as the stop sector is not measured at the LHC. 
Therefore the study is restricted to the parameters of the 
first two generations. 

The bottom-up evolution of the scalar sector for the 
first two generations in DS7 is shown in Figure [S] The 
unification is qualitatively observed at about 10^^ GeV as 
expected. While the slepton parameters are measured pre- 
cisely at the electroweak scale, the coupling of the RGBs 
leads to a quick degradation of their precision as function 
of the scale. It is obvious that the scalars will not be able 
to improve the determination of the unification scale and 
therefore will not improve the separation of DSl and DS7. 



5.2.1 Unification scale and unified parameters 

Given the observation of unification in the gaugino and 
scalar sector, the unification scale and the unified param- 
eter can be determined at the LHC. Here the study is 
restricted to the true solution (DS7) without loss of gen- 
erality, the solution of DSl leads to a similar precision. 
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Fig. 5. Evolution of the first and second generation scalar mass 
parameters for the true solution (DS7) at the LHC : bottom-up 
evolution of the MSSM. 



Name Unified Unification Scale Parameter at 

Parameter [log(Q/GeV)] L7 -lO" GeV 



mi/2 
mo 



251.9±5.9 
98.5±10.5 



16.23±0.29 
16.5±0.6 



252.3±3.2 
100.8±4.9 



Table 6. Measurement of grand unification with LHC mea- 
surements (DS7). All masses are in GeV. 



The results for the gaugino mass parameter mi/2 as 
well as the scalar mass parameter mo are shown in Table[S] 
For the trilinear couplings the calculation is not useful as 
only one parameter is free (At) and the other two are fixed 
at the electroweak scale. 

The most precise determination of the unification scale 
is obtained in the gaugino sector with a measurement of 
(1.7±1.1) -10^^ GeV. At the unification scale mi/2 is mea- 
sured to about 2% and is in agreement with the nominal 
value of SPSla (250 GeV). Fixing the scale reduces the er- 
ror on the common mass parameter by almost a factor 2. 

The common scalar parameter mg is determined with 
a precision of about 10% in agreement with the nominal 
value of SPSla of 100 GeV. As the scale is measured more 
precisely in the gaugino sector, combining the two sectors 
will not provide an improvement. Alternatively one can 
determine mo at this fixed scale: mo is measured to be 
101 GeV with an error of 5 GeV, thus the error is re- 
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duced by a factor 2, not including the error on the scale 
determination. 

Thus once the ambiguous solutions for the LHC are 
discarded, the common scalar and gaugino mass parame- 
ters can be reconstructed in a bottom-up approach with 
a precision of 10% and 2% respectively. The precision is 
improved to and 5% and 1% respectively if the unification 
scale is fixed. The fixed scale results can be compared 
to the mSUGRA parameter determination reported in 
Ref. [30] where a precision of 2% was reached on the scalar 
mass and roughly 1% on the on the gaugino mass. Thus 
in the gaugino sector the results of top-down (mSUGRA) 
and bottom-up (MSSM) agree well, whereas in the scalar 
sector the determination bottom-up is less precise. 



5.2.2 Effect of threshold corrections at the high scale 

A further complication arises from threshold corrections 
of the unknown embedding theory. The size and the sign 
of such corrections are model-dependent, but in typical 
SU(5) models these primarily affect [641165] as{McuT) 
(and Mgut to some extent). This could compensate for 
the observed mismatch, at two-loop RGE in MSSM, in 

as{MGUT) - a2{MGUT) ^ (a2(MG;7T) = ai{MGUT) = 
gf/{4:Tr)). For a typical minimal SUGRA input, the latter 
mismatch is a few percent and negative (about -3% in par- 
ticular for SPSla). In addition, intrinsic corrections to the 
gaugino masses (i.e. corrections to Eq. [5]) have been eval- 
uated to be a few percent in a minimal SU(5) mo del [551. 
i.e. roughly of the order of two-loop MSSM corrections, 
though the former can be much larger in non-minimal 
GUT models, e.g. with large representations of heavy chi- 
ral multiplets. Since specific GUT model corrections are 
anyway not implemented at present in the spectrum calcu- 
lators, for simplicity a positive shift in M3 correlated with 
the as one (i.e. preserving Eq. [5]at one-loop) is assumed. 
The effect is thus approximated by shifting the measured 
parameter M3 by 3% for illustration, while possible model- 
dependent effects on other parameters are neglected. 

The parameter mi/2, including the threshold correc- 
tions, is shifted by 3.5 GeV and the unification scale by 
0.07, corresponding to a shift of 0.3 • 10^^ GeV. The abso- 
lute values of the shifts have to be compared to the error 
of the determination of the common mass (5.9 GeV) and 
the scale (0.3). The shift corresponds to a deviation of 
less than half a standard deviation for the mass, thus the 
threshold effects will not play a large role at the LHC, 
given the expected precision. 

In addition to the study of the central values, it is also 
interesting to address the question whether the thresh- 
old effects could lead to the conclusion that DSl unifies 
and DS7 does not. The percentage of the unified param- 
eter sets at the best scale is a good indicator. Including 
the threshold corrections, in DS7 the percentage drops to 
87.4% (from 95.4%). In DSl, the unification percentage in 
the gaugino sector is 3% (from 38%). Thus the threshold 
corrections applied to DSl will actually increase the differ- 
ence between the true and the wrong solution. However, 
if the sign of the M3 shift is opposite (i.e. if the specific 



GUT model is such that those corrections are larger and 
essentially negative) this conclusion would change. 



5.2.3 Evolution with shifted data 

The studies described so far all dealt with datasets which 
are smeared, but centered around the true central value. 
An additional complication will arise with real data as the 
measured value of the parameters will be shifted, within 
the theoretical and experimental errors, from the true cen- 
tral value. In this case it is still possible to use the toy 
experiments but they are performed around the shifted 
values. 

It is also necessary to verify that the theoretical er- 
rors used in the study cover at least the difference of the 
predictions from spectrum generators which have similar 
precision. The dataset (SPSla) calculated by SUSPECT is 
used, but the MSSM parameters are determined by using 
SoftSUSY, i.e., SoftSUSY is used to predict the spectrum 
and evolve the parameters to the high scale. The SPSla 
dataset from SUSPECT corresponds to a shifted dataset 
for SoftSUSY. 

The common gaugino mass parameter mi/2 is deter- 
mined to be 252.7 ± 6.4 GeV at log(g/GeV) = 16.2 ± 0.3 
with SoftSUSY. The difference with respect to the de- 
termination using only SUSPECT is less than about one 
standard deviation (Table[B]). In the scalar sector the com- 
mon mass is determined to be 92 GeV with an error of 
10 GeV. The results, both the central values and errors, 
at the EW scale as well as the GUT scale, are in excel- 
lent agreement with SUSPECT for the gaugino and scalar 
sector at the LHC. 

It is also interesting to note that the percentage of 
toy sets compatible with grand unification in the gaugino 
sector is essentially unchanged: 95% for DS7, 35% for DSl. 

The results show that as required by the definition of 
the size of the theoretical error at least the difference be- 
tween different spectrum calculators is covered. While dif- 
ferences will be observed, depending on which calculator 
is used, the difference is small with respect to the expected 
error on the parameters and the unification scale. 



5.2.4 High scale MSSM 

To study the difference between bottom-up and top-down 
evolution, 5000 toy experiments are used to determine the 
parameters of the high scale MSSM. The results are com- 
pared to the parameters of the MSSM determination after 
their evolution to the high scale. Parameter sets are re- 
moved as soon as they became tachyonic. Thus they con- 
tribute to the RMS below the scale where they become 
tachyonic, but not above. 

The precision of the parameters is comparable in the 
gaugino sector at the level of 20-30%. In the scalar sec- 
tor the precision of the high scale MSSM parameters is 
reduced by a factor of more than 5. The difference with 
respect to the bottom-up evolution (Figure [5]) is obvious 
in Figure [S] 
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parameters compared to the high scale MSSM where such 
parameter sets are excluded by construction. 



5.3 Evolution of the parameters from LHC+ILC 
observables 

In Table \5\ the result of the determination of the MSSM 
parameters is shown in the third column for the combina- 
tion of LHC and ILC. With the exception of the tri-linear 
couplings, the parameters are measured with excellent pre- 
cision at the clcctroweak scale. Additionally the eight-fold 
ambiguity left by the LHC data alone is solved by the ILC. 



Fig. 6. Evolution of the first and second generation scalar mass 
parameters for the true solution (DS7) at the LHC : top-down 
evolution for the high scale MSSM. 



The apparent contradiction between the naive expec- 
tation that top-down should be equal to bottom-up and 
the large differences that are observed, can be traced to 
several factors: in contrast to the top-down evolution where 
all parameters are defined at a single high scale, in the 
bottom-up case initial conditions for the RGE are more 
sensitive to the (many) different physical scales. All thresh- 
old corrections are calculated at present in the one-loop 
approximation for most sparticles. This induces an in- 
creased sensitivity to the intrinsic errors. 

In the stau sector only one measurement is available 
at the LHC. This measurement can be used to determine 
Mfj^. In the MSSM the parameter Mf^ decouples and can 
take almost any value, even 1 TeV. This is reflected in 
the large error in Table [S] However, in the high scale 
MSSM, the top-down running introduces interdependen- 
cies. Keeping all parameters at their nominal SPSla value 
and moving only M^^ to 200 GeV is not possible: the 
selectron and smuon masses are changed by 5 GeV. As 
the lepton-lepton edge, which depends on these slepton 
masses, is measured precisely at the LHC, such a large 
change is not compatible with the observables. Thus the 
high scale MSSM restricts the nominally available param- 
eter space. 

All scalar parameters are correlated non-linearly through 
Tr[yTO^] entering all RGE equations. In a bottom-up evo- 
lution, large errors on the MSSM sfermion parameters 
induce a departure from zero in the boundary (initial) 
value of Tr[yTO^]Since the dependence is quadratic, this 
is strongly amplified in the evolution. In the top-down evo- 
lution as the sfermion parameters are interdependent due 
to the RGBs and as a consequence all parameters are well 
measured, the departure of the Tr[ym^] from its SPSla 
value is less pronounced and its effect smaller. 

While the tachyonic parameter sets are removed and 
play no role at the high scale, they contribute to the er- 
ror on the parameters at all scales below the GUT scale. 
These parameter sets are relatively far away from their 
nominal SPSla value and thus lead to larger RMS for the 
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[log(0/GeV)] 


2.33 -lO^^ GeV 


mi/2 


249.5±1.8 


16.37±0.05 


249.6±1.5 


l/2Gcn 

mo 


98.2±10.7 


16.5±0.7 


100.4±2.5 


-„3Gcn 

mo 


117.1±27 


15.4±L3 


103.1±25 


mo 


105.3±9.1 


15.9±0.6 


99.4±2.0 


Ao 


-164±182 


14.8±4.5 


-133.8±207 



Table 7. The results for the measurement of the common pa- 
rameters and unification scale with LHC+ILC measurements 
in the bottom-up approach are shown. All masses are in GeV. 
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Fig. 7. The evolution of the gaugino masses is shown for the 
combined results from LHC-I-ILC. 



The results of the evolution of the parameters mea- 
sured at the electroweak scale as well as the unification 
are shown in Table [71 The evolution of the three gaugino 
mass parameters is shown in Figure [T] 

The common gaugino parameter mi/2 is determined 
with a precision of 1.8 GeV at a grand unification scale 
of (2.33 ± 0.28) ■ IQi^ GeV in agreement with the SPSla 
parameter set definition. Note that while the logarithm 
of the unification scale is determined with a precision of 
0.3%, due to the proper error propagation, the scale in 



Claire Adam et al.: Measuring Unification 



13 



GeV is only determined with 10% precision. The error 
with respect to LHC alone is reduced by a factor 3. 

Fixing the unification scale to the central value de- 
termined by the gauginos, the error on mi/2 would be 
reduced by about 0.2 GeV. It is also instructive to ana- 
lyze the individual contributions of the three parameters 
to the determination of the unified mass parameter. At 
the unification scale Mi, M2 and M3 are measured with a 
precision of 1.6 GeV, 2.0 GeV and 3.3 GeV respectively. 
Thus the error on m^ /2 is essentially equal to the preci- 
sion of Ml. The naive combination of the three parameter 
errors, i.e., ignoring correlations, would lead to an error 
of 1.1 GeV, almost 30% better than that obtained. While 
M3 is uncorrelated with the other two parameters. Mi and 
M2 are almost 100% positively correlated. Thus combin- 
ing the latter two will not increase the precision. M3 on 
the other hand, while not correlated, is less precisely mea- 
sured and therefore the combined error is decreased only 
by a small amount in the combination as the error (cr) on 
the combination of two uncorrelated measurements (cti, 
0-2) reads: 



large negative correlations among the parameters which 
reduce the total error. 




1/Jl/aj + l/al 



(16) 



Fig. 9. The evolution of the scalar parameters of the third 
generation is shown for the combined results from LHC-f ILC 
as function of the logarithm of the scale. 
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Fig. 8. The evolution of the scalar masses of the first genera- 
tion is is shown for the combined results from LHC-I-ILC. 



Due to the increased precision of LHC-I-ILC, the uni- 
fication of the scalar mass parameters can be separated 
into the first two generations and the 3rd generation. The 
evolution of the scalar masses of the first generation is 
shown in Figure [S] As shown in Table [71 iuq "^^ can be 
determined with a precision of about 10% in agreement 
with the nominal value (100 GeV) of SPSla. The deter- 
mination of the logarithm of unification scale is less pre- 
cise than the precision in the gaugino sector by an order 
of magnitude. Using the unification scale defined by the 
gaugino measurement, the error is reduced to 2.5 GeV. 
It is interesting to note that in the scalar case the naive 
combination of the parameters neglecting the covariance 



matrix would lead to an error on 



l/2Ge 



of 8.2 GeV, 



thus greater than the correct value. This is due to the 



For the 3rd generation the unification procedure leads 
to a less precise determination (27 GeV) of the common 
scalar parameter as can be immediately inferred from Fig- 
ure [9] as well as Table [T] The improvement by fixing the 
unification is rather small. The larger errors with respect 
to the first two generations are due to two sources. The 
parameters of the third generations are less precisely mea- 
sured than those of the first two generations. Additionally, 
as discussed in Section [2Tl the sfermion mass terms in the 
RGE have a stronger inter-dependence. Nevertheless the 
reconstructed unification scale is in agreement with the 
SPSla parameter set. 

If instead of separating the first and second generation 
from the third generation (and the two Higgs parame- 
ters), all parameters relating to scalars are combined, the 
error on the determination of the unification scale as well 
the parameter is decreased as shown in the second to last 
row of Table [T] Using the central value of the unification 
scale determined by the gaugino sector, the error on mg 
is reduced to 2 GeV. This error is smaller than the naive 
combination ignoring correlations by a factor 3, showing 
the necessity of a proper treatment of the errors and cor- 
relations. 

It is interesting to note that the error here is identical 
to the error from the LHC alone. This might seem surpris- 
ing at first sight as the slepton sector is measured experi- 
mentally an order of magnitude more precisely at the ILC 
than at the LHC. The reason for this (superficial) lack of 
impact lies again in the structure of the RGEs. The error 
on the sleptons, as shown in Figure [51 increases strongly 
as function of the scale due to the coupling with the less 
precisely determined squark sector. Thus at the unifica- 
tion scale, the slepton precision of the ILC is diluted by 
the LHC squark precision. Additionally while for the com- 
bination of LHC and ILC no MSSM parameters are fixed, 
for the LHC the trilinear couplings of the sbottoms and 
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staus are fixed. Fixing these leads to contraction of the al- 
lowed parameter space and therefore an artificial decrease 
of the scalar mass error at the LHC. 
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Fig. 10. (left) The evolution of the trilinear parameters of 
the third generation is shown for the combined results from 
LHC+ILC. (right) The evolution of the squares of the Higgs 
sector parameters Mhi and Mh2 is shown for LHC+ILC. 



Of the trilinear couplings, only At can be measured 
with good precision with the mass measurements from 
the ILC and the edges from the LHC. Further measure- 
ments will be necessary to constrain these parameters. For 
completeness sake the results are shown in the last line of 
Table [7] and in Figure [TUl (left). Due to the large errors on 
the parameters determined at the electroweak scale, the 
parameters are compatible with unification at all scales. 
Fixing the unification scale to the one determined by the 
gaugino sector, the error on Aq is 200 GeV. This is a slight 
improvement compared to the error of 211 GeV on At at 
the unification scale. Only A,- contributes to the reduc- 
tion of the error as A{, has a larger error and is therefore 
irrelevant for the combination. 

The evolution of the squares of the Higgs sector param- 
eters Mhi and Mhs , replacing the electroweak parameters 
ruA and /j, is shown in Figure [TU] (right). The square of 
Mh2 is negative at the electroweak scale, as required by 
electroweak symmetry breaking, and reasonably well de- 
termined. As the evolution proceeds to higher scales, the 
error increases significantly, somewhat faster for Mhi . At 
the high scale the parameters are compatible with the true 
value of (100 GeV)^, but also with zero. Thus there is no 
significant contribution on the determination of nip from 
these parameters. 

The fixed scale results of the common parameters at 
the high scale can be compared to the mSUGRA param- 
eter determination reported in Ref. [5^. While for the 
gauginos the error of the bottom-up determination com- 
pared to mSUGRA is about 2.5 times larger, the common 
scalar mass is determined with a precision 5 times better 
in mSUGRA. The largest difference is observed for the 
tri-linear coupling where the mSUGRA determination is 
more precise than the bottom- up one by a factor of 20. 
This shows that bottom-up and top-down do not lead to 
the same results. 

The last question to be addressed is to determine the 
probability with which grand unification will be measured 
in the four measurements. At the unification scale 95% of 
the toy experiments show unification in the gaugino sec- 
tor in agreement with the definition of the x^ cut. For 



the scalar parameters at the unification scale defined by 
the gauginos, the most precise measurement of this scale, 
90% unify. For the trilinear couplings 93% of the toy ex- 
periments are compatible with grand unification. 

The additional observables added by the ILC to the 
LHC dataset indeed increase the precision of the determi- 
nation of the couplings in the gaugino sector by a large 
factor. Additionally different unification hypotheses can 
be tested (full scalar unification, separate unification for 
the light and heavy generations). 



5.3.1 Effect of threshold corrections at the high scale 

To study the effect of threshold corrections the measured 
value of M3 is shifted by 3% as in Section [5. 2. 2 1 The com- 
mon gaugino mass is then determined to be 251.7 GeV 
and the scale is shifted by 0.05. The mass shift at the 
LHC alone is less than about half of a sigma. Due to the 
higher precision of the combination LHC-I-ILC, the shifts, 
while similar in absolute numbers, is now of the order of 
a sigma. 

The effect on the percentage of unifying parameter sets 
is also much larger: it decreases from 95% to 77%. For 
comparison, the effect is only half as large for the LHC 
alone. The increased precision added by the ILC means 
that the effect of threshold effects become more important. 



5.3.2 Evolution with shifted data 



As in Section 15.2.31 to illustrate the effect of a shifted 
dataset, the central values from SUSPECT are used, but 
the predictions as well as the evolution bottom-up are per- 
formed by SoftSUSY. 

In the gaugino sector the unification scale is deter- 
mined with the same precision as before. The central val- 
ues differ by less than 0.01, corresponding to one fifth of 
a sigma. mi/2 is determined to be 249.5 GeV, in excel- 
lent agreement with the determination using SUSPECT. 
The error on mi/2 is comparable at 1.6 GeV. 94% of the 
parameter sets unify. 

In the scalar sector mo is shifted by 1.6 GeV closer to 
the nominal value of SPSla. The error from the determi- 
nation of 9 GeV is comparable to that from SUSPECT. 
Thus the shift is less than one fifth of a sigma. 

The study shows that the theoretical errors fulfill the 
requirement that they cover at least the difference between 
different calculations of the spectrum and the RGE run- 
ning. The results are in excellent agreement, showing that 
the analysis is robust. 



5.3.3 High scale MSSM 

As shown in Figure[l] (right), the fraction of "non-tachyonic" 
datasets is much higher in the LHC-I-ILC case. But the 
comparison of the errors on the MSSM parameters evolved 
to the high scale with the high scale MSSM parameters 
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shows again the difference between a top-down and a bottom- 
up approach. The difference between the errors on the pa- 
rameter determination at the unification scale is reduced 
due to the increased precision added by the ILC measure- 
ments. In the gaugino sector these are now of the order of 
10-20%. In the scalar sector the differences remain much 
larger, a factor 5 in the stau sector for example. The re- 
sults of the top-down evolution shown in Figure [TT] can 
easily be compared with the bottom-up evolution shown 
in Figure [To] (right) and Figure [9] 
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Fig. 11. The top-down evolution of high-scale MSSM param- 
eters at the LHC-I-ILC is shown for : (left) third generation 
scalar mass parameters, (right) squares of the Higgs sector pa- 
rameters Mhi and Mhj- 



also far away from their nominal SPSla values at the high 
scale. 

The x^ of both parameter determinations is very good, 
less than the degrees of freedom with a slightly smaller 
value for the MSSM bottom-up resuh {Ay^ = 3/18d.o.f.). 
The result of the evolution of the two high scale MSSM 
determinations is shown in Figure [12] for the third genera- 
tion scalar parameters and Mhs • All parameters with the 
exception of M^ and Mqa^^ (and the tri-linear couplings) 
converge from extremely different values at the high scale 
to the same value at the EW scale. This effect can be un- 
derstood from the fact that at the leading one-loop RGE, 
both Mjj^ and Mqs^ evolution depend strongly on Mhs, 
which exhibits a strong variation around the EW scale as 
mentioned above. In contrast, other relevant scalar masses 
do not depend (at one-loop RGE level) on Mhj • Moreover 
the other scalar parameters RGE depend very little on 
Mjj^ and Mqs^, the effect is suppressed by the bottom 
Yukawa coupling, which is quite small for SPSla. Indeed, 



the evolution of the combination 2 • M-g 



M? , which 



eliminates the dominant dependence on Mfja at one-loop 
RGE, has essentially the same form for the two solutions 
as expected. 



The strong sensitivity to the Higgs mass term Mhs as 
illustrated in Figure 1111 is not surprising since in most 
MSSM scenario its evolution, mainly driven by the top 
Yukawa coupling, drastically accelerates near the EW scale, 
where generally M^ changes sign. In other words the 
slope of its beta function becomes large around the EW 
breaking scale, which is taken as the initial scale in a 
bottom-up evolution. Therefore in the bottom-up evolu- 
tion even a small error in the initial low scale Mh2 value 
can induce a very large difference at the high scale. In 
contrast, in a top-down evolution near the GUT scale the 
initial slope of the beta function for Mjja is moderate, and 
the final value at the EW scale is less dependent on ini- 
tial high scale boundary conditions, illustrating a typical 
focusing behavior. 

To illustrate the impact of this behavior on the scalar 
mass determination, a bottom-up parameter set is selected 
for which Mgj^, after evolution to the high scale, is about 
60 GeV, i.e., far away from the SPSla nominal values. 
The choice of Mgj^ to select the dataset/parameter set is 
motivated by the fact that it is well measured and Mg^ is 
the lightest scalar mass. This parameter is therefore more 
likely to evolve to tachyon values in a bottom-up evolu- 
tion. A different choice of the scalar parameter departing 
from its SPSla value would also have been possible. 

Two high scale MSSM parameter determinations are 
performed with this dataset: one starting from the nom- 
inal SPSla values and a second one starting from the 
MSSM results after evolution to the high scale. In the 
first case, the high scale MSSM scalar masses are deter- 
mined at values within 2 GeV of the SPSla nominal values 
for the slepton masses. The best-fit result of the second 
determination is compatible with the result of the bottom- 
up study, where apart from Msj, other scalar masses are 
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Fig. 12. The top-down evolution of high-scale MSSM param- 
eters at the LHC-I-ILC is shown for ; (left) third generation 
scalar mass parameters, (right) squares of the Higgs sector pa- 
rameters Mhi and Mh2 • The two lines for each parameter cor- 
respond to the two essentially degenerate solution for the same 
dataset. One solution is close to the SPSla values at the high 
scale, the other one far away. At the EW scale the mass pa- 
rameters with the exception of Mj and MqSj^ converge to the 
same value. 



Given the smallness of Ay^ ^ the two solutions are de 
facto degenerate. It is interesting to ask why the bottom- 
up solution has found an ever so slightly better solution 
than the standard high scale study. While at the EW scale 
the difference of Mj^ and Mqs^^ for the two solutions is of 
the order of 20 GeV, i.e., close by, at the high scale the 
difference is of the order of 150 GeV. Thus at the EW 
scale the sampling of the parameter space is much easier: 
all mass parameters but two are the same and only a small 
excursion of 20 GeV in Mj^ and Mq3^ is needed to find and 
differentiate between the two solutions. At the high scale 
however, all scalar parameters of the two solutions are far 
apart: 60GeV for Mg^, , 150 GeV for Mj^ etc. Therefore the 
parameter determination will find easily the solution close 
to the values of SPSla, while it is more difficult to find 
the second solution as a much larger parameter space has 
to be sampled. But even in this case, by construction, the 
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top-down extrapolation will miss the contribution of the 
tachyonic parameter sets from the bottom-up approach. 

The high scale MSSM top-down results presented in 
this study agree with the bottom-up results presented by 
other groups |261I67) . In Refs. p6 l [g7 ) additional ILC ob- 
servables are used (polarized cross sections) and the theory 
errors are set to zero. Additionally in Ref. [2S] the trilin- 
ear parameters A^ and A^ are required to be at the high 
scale compatible with At within 2 sigma. This additional 
requirement leads mechanically to a reduction of the error 
on the corresponding mass parameters. 
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The discovery of supersymmetry at the LHC will lead to a 
wealth of signatures which can be exploited to determine 
many MSSM parameters. In parameter regions similar to 
the SPSla parameter point, they can be determined at 
the LHC up to an at least 8-fold ambiguity in the gaug- 
ino sector. Although a part of those ambiguities may be 
resolved, e.g., by a complementary study of the MSSM 
contributions to the dark matter relic density, a full res- 
olution will likely require a complete observation of the 
sparticle spectrum at the ILC. 

Starting from the electroweak scale, we can test the 
unification of different supersymmetry-breaking parame- 
ters. While remaining ambiguities make it impossible to 
measure unification at the LHC, it will nonetheless be 
possible to classify solutions into the ones compatible and 
the ones not compatible with unification. In the case of 
an ambiguous solution (DSl) which differs from the true 
solution only by the sign of /x, the differentiation will be 
difficult as about 38% of the parameter sets corresponding 
to this wrong solution nevertheless unify. 

This way, at the LHC the unified gaugino mass pa- 
rameter can be measured bottom-up to about 2% and the 
logarithm of the unification scale to 1.7%. Adding the ILC 
data improves the determination of the mass by more than 
a factor 3 and the unification scale by almost one order of 
magnitude. In the scalar sector the errors are generically 
larger at the level of 10%. The errors on the trilinear cou- 
plings are too large to be used for a determination of the 
unification scale. 

The robustness of our results we have confirmed by 
comparing two different renormalization group tools: SUS- 
PECT and SoftSUSY. The parameter determination as 
well the evolution are in good agreement within the er- 
rors. Threshold corrections at the high scale were studied 
for a particular model, motivated by SU(5) grand unifica- 
tion. The percentage of parameter sets unifying is affected 
more strongly including the ILC observables providing an 
increased sensitivity. 

Finally, our study show that a proper bottom-up ap- 
proach will clearly lead to different results from simply 
determining the parameters of the high scale MSSM (or 
mSUGRA). In addition to resolving the ambiguities at the 
LHC, the ILC plays a strong role in the stabilization of 
the validity of the parameter sets as function of the scale. 
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